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Abstract 

An amplitude analysis of the final state structure in the B^^ — )• J /iIjK^ K~ decay 
mode is performed using 1.0 fb~^ of data collected by the LHCb experiment in 7 TeV 
center-of-mass energy pp collisions produced by the LHC. A modified Dalitz plot 
analysis of the final state is performed using both the invariant mass spectra and the 
decay angular distributions. Resonant structures are observed in the K^K~ mass 
spectrum as well as a significant non-resonant S-wave contribution over the entire 
K^K^ mass range. The largest resonant component is the 0(1020), accompanied 
by /o(980), /2(1525), and four additional resonances. The overall branching fraction 
is measured to be B{Bl J/i^K+K-) = (7.70 ± 0.08 ± 0.39 ± 0.60) x 10-^ where 
the first uncertainty is statistical, the second systematic, and the third due to the 
ratio of the number of to B~ mesons produced. The mass and width of the 
/^(1525) are measured to be 1522.2 ± 2. 8l|;j MeV and 84 ±6^^^ MeV, respectively 
The final state fractions of the other resonant states are also reported. 
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1 Introduction 

The study of decays to J/%l)h^h~ , where h is either a pion or kaon, has been used to 
measure mixing-induced CP violation in 5° decays hjiTlFl In order to best exploit these 
decays a better understanding of the final state composition is necessary. This study has 
been reported for the 5° — )■ J/TpTi^Ti~ channel js]. Here we perform a similar analysis for 
5° J/^jjK^K~. While a large 0(1020) contribution is well known jo] and the /2(1525) 



component has been recently observed [10] and confirmed 11 , other components have not 
heretofore been identified including the source of S-wave contributions [12]. The tree-level 
Feynman diagram for the process is shown in Fig. [T| 




Figure 1: Leading order diagram for — J/tpK^K . 

In this paper the J/ ipK^ and K^K~ mass spectra and decay angular distributions are 
used to study resonant and non-resonant structures. This differs from a classical "Dalitz 
plot" analysis [13j since the J/tfj meson has spin-1, and its three helicity amplitudes must 
be considered. 



2 Data sample and detector 

The event sample is obtained using 1.0 fb^^ of integrated luminosity collected with the 



LHCb detector 14 using pp collisions at a center-of-mass energy of 7 TeV. The detec- 
tor is a single- arm forward spectrometer covering the pseudorapidity range 2 < 77 < 5, 
designed for the study of particles containing 6 or c quarks. Components include a high 
precision tracking system consisting of a silicon-strip vertex detector surrounding the pp 
interaction region, a large-area silicon-strip detector located upstream of a dipole magnet 
with a bending power of about 4Tm, and three stations of silicon-strip detectors and 
straw drift-tubes placed downstream. The combined tracking system has momentum]^ 
resolution Ap/p that varies from 0.4% at 5GeV to 0.6% at 100 GeV. The impact param- 
eter (IP) is defined as the minimum distance of approach of the track with respect to the 
primary vertex. For tracks with large transverse momentum with respect to the proton 
beam direction, the IP resolution is approximately 20 \\m.. Charged hadrons are identified 
using two ring-imaging Cherenkov detectors. Photon, electron and hadron candidates are 



'''Mention of a particular mode implies use of its charge conjugate throughout this paper. 
^We 'work in units 'where c = 1. 
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identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. 

The trigger 15 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage that applies a full event reconstruction. 
Events selected for this analysis are triggered by a J/V' — )■ /^^A'-" decay, where the J /iIj is 
required at the software level to be consistent with coming from the decay of a meson 
by use either of IP requirements or detachment of the J/ip from the primary vertex. 
Monte Carlo simulations are performed using Pythia |16| with the specific tuning given 
in Ref. 17 , and the LHCb detector description based on Geant4 (Ts] described in 
Ref. IT9]. Decays of B mesons are based on EvtGen 120|. 



3 Signal selection and backgrounds 

We select 5^ — t- J /ipK^ K~ candidates trying to simultaneously maximize the signal 
yield and reduce the background. Candidate J/ip ^ /^^A*" decays are combined with 
a pair of kaon candidates of opposite charge, and then requiring that all four tracks 
are consistent with coming from a common decay point. To be considered a J/ip ^ 
fj,'^fj,~ candidate, particles identified as muons of opposite charge are required to have 
transverse momentum, p^, greater than 500 MeV, and form a vertex with fit P^i' 
number of degrees of freedom (ndf) less than 11. These requirements give rise to a large 
J/ip signal over a small background [21] . Only candidates with a dimuon invariant mass 
between —48 MeV to +43 MeV relative to the observed J/ip mass peak are selected. The 
asymmetric requirement is due to final-state electromagnetic radiation. The two muons 
are subsequently kinematically constrained to the known J/ip mass [o]. 

Our ring-imaging Cherenkov system allows for the possibility of positively identifying 
kaon candidates. Charged tracks produce Cherenkov photons whose emission angles are 
compared with those expected for electrons, pions, kaons or protons, and a likelihood for 
each species is then computed. To identify a particular species, the difference between the 
logarithm of the likelihoods for two particle hypotheses (DLL) is computed. There are 
two criteria used: loose corresponds to DLL{K — 71) > 0, while tight has DLL(K — vr) > 10 
and DLL{K — p) > —3. Unless stated otherwise, we require the tight criterion for kaon 
selection. 

We select candidate K^K~ combinations if each particle is inconsistent with having 
been produced at the primary vertex. For this test we require that the formed by using 
the hypothesis that the IP is zero be greater than 9 for each track. Furthermore, each 
kaon must have > 250 MeV and the scalar sum of the px of the kaon candidates must 
be greater than 900 MeV. To select 5° candidates we further require that the two kaon 
candidates form a vertex with < 10, and that they form a candidate vertex with 
the J /if) where the vertex fit x^/^df < 5. We require that this vertex be more than 
1.5 mm from the primary vertex, and the angle between the 5° momentum vector and 
the vector from the primary vertex to the 5° vertex must be less than 11.8 mrad. 
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The candidate invariant mass distribution is shown in Fig. |2] The vertical hnes 
indicate the signal and sideband regions, where the signal region extends to ±20 MeV 
around the nominal mass [o] and the sidebands extend from 35 MeV to 60 MeV on 
either side of the peak. The small peak near 5280 MeV results from B^ decays, and will 
be subject to future investigation. 
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Figure 2: Invariant mass spectrum of J /'ipK^K combinations. The vertical lines indicate the 
signal (black-dotted) and sideband (red-dashed) regions. 

The background consist of combinations of tracks, which have a smooth mass shape 
through the J/ipK^K' region, and peaking contributions caused by the reflection of 
specific decay modes where a pion is misidentified as a kaon. The reflection background 
that arises from the decay B^ J /ipK'TT^ , where the vr"*" is misidentified as a , is 
determined from the number of B^ candidates in the control region 25 — 200 MeV above 
the B^g mass peak. 

For each of the candidates in the J/ipK~^K~ control region, we reassign each of the 
two kaons in turn to the pion mass hypothesis. The resulting J/ipKn invariant mass 
distribution is shown in Fig. [3] The peak at the B^ mass has 906 ± 51 candidates, 
determined by fitting the data to a Gaussian function for the signal, and a polynomial 
function for the background. From these events we estimate the number in the B^ signal 
region, based on a simulation of the shape of the reflected distribution as a function of 
J/^K-K+ mass. Using simulated 5° J/?/^F*°(892) and B^ J/V'F*(1430) samples, 
we calculate 309 ± 17 reflection candidates within ±20 MeV of the 5° peak. This number 
is used as a constraint in the mass fit described below. 

To determine the number of 5° signal candidates we perform a fit to the candi- 
date J /il>K^K~ invariant mass spectrum shown in Fig. |4j The fit function is the sum 
of the -Bg signal component, combinatorial background, and the contribution from the 
B^ — J/il)K~Ti^ reflections. The signal is modeled by a double- Gaussian function with 
a common mean. The combinatorial background is described by a linear function. The 
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Figure 3: Invariant mass distribution for J /ifjK^ K~ candidates 25 — 200 MeV above the 
mass, reinterpreted as — )■ J/'ipK^n^ events. The fit is to a signal Gaussian whose mass and 
width are aUowed to vary as well as the polynomial background. 




m(J/ilJ K"K") (MeV) 



Figure 4: Fit to the invariant mass spectrum of J /tpK'^ K~ combinations. The dotted (black) 
line is the combinatorial background, the dashed (red) shape shows the misidentified — )• 
J/^/'i^~7r+decays, and the solid (blue) curve shows the total. The vertical dashed lines indicate 
the signal region. 



reflection background is constrained as described above. The mass fit gives 19,195±150 
signal together with 894 ± 24 combinatorial background candidates within ±20 MeV of 
the 5° mass peak. 

We use the decay B" — )■ J/ipK' as the normahzation channel for branching fraction 
determinations. The selection criteria are similar to those used for J/ip K~^K~ , except for 
particle identification as here a loose kaon identification criterion is used. Figure [5] shows 
the J/ipK~ mass distribution. The signal is fit with a double- Gaussian function and a 
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linear function is used to fit the combinatorial background. There are 342,786±66f signal 
and 10,195±134 background candidates within ±20 MeV of the B~ peak. 




5250 5300 

m(J/i|j K") (MeV) 

Figure 5: Fit to the invariant mass spectrum of J/'ipK~ candidates. The dotted Une shows the 
combinatorial background and the solid (blue) curve is the total. 



4 Analysis formalism 

One of the goals of this analysis is to determine the intermediate states in — > 
J /'il)K^K~ decay within the context of an isobar model (22p3] , where we sum the resonant 
and non-resonant components testing if they explain the invariant mass squared and angu- 
lar distributions. We also determine the absolute branching fractions of — > J/V' 0(1020) 
and — t- J/ip f 2(1525) final states and the mass and width of the /2(1525) resonance. 
Another important goal is to understand the S-wave content in the 0(1020) mass region. 

Four variables completely describe the decay of 5° — )• J/ipK^K~ with J/ip ^ fi'^fj,~. 
Two are the invariant mass squared of J/ipK^, = m?{J /ipK^), and the invariant 
mass squared of K~^K~, S23 = m'^{K^K'). The other two are the J/ijj helicity angle, 
6j/^, which is the angle of the fi~^ in the J/ip rest frame with respect to the J/ip direction 
in the B^ rest frame, and the angle between the J/ip and K^K~ decay planes, Xi the 
5^ rest frame. To simplify the probability density function (PDF), we analyze the decay 
process after integrating over the angular variable x? which eliminates several interference 
terms. 

4.1 The model for Jj^K^K' 

In order to perform an amplitude analysis a PDF must be constructed that models cor- 
rectly the dynamical and kinematic properties of the decay. The PDF is separated into 
two components, one describing signal, S", and the other background, B. The overall PDF 
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given by the sum is 

F{Si2,S23,9j/^) = '^£{Sl2,S23,0j/^)S{Si2,S23,9j/,^) (1) 

+ B{si2,S23,dj/-4,), 

where e is the detection efficiency. The background is described by the sum of combina- 
torial background, C, and reflection, R, functions 

B{Si2, S23, 0J/ip) = -^P^C(Si2, S23, Oj/^) + ■j^R{Si2, S23, dj/i;), (2) 

where /com and fren are the fractions of the combinatorial background and reflection, 
respectively, in the fitted region. The fractions /com and /res obtained from the mass fit 
are fixed for the subsequent analysis. 

The normalization factors are given by 



A4ig = J £{Sl2,S23,0j/^)S{Si2,S23,0j/i,)dSi2dS23dcOs6j/^, 
A/'com = j C{Si2,S23,0j/^) dSi2dS23dcOs9j/^, (3) 
A/'rcfi = j R{Sl2,S23,0j/^) dSi2dS23dcOs9j/^. 



This formalism similar to that used by Belle in their analysis of — )■ K^n^Xci |24j, and 
later used by LHCb for the analysis of B^ — )■ J/ ^/iTr+Tr" [s] . 

The invariant mass squared of J/i/jK~^ versus K^K~ is shown in Fig. [g] for 5^ — > 
J/ipK^K^ candidates. No structure is seen in m?{J/il) K^). There are however visible 
horizontal bands in the K^K~ mass squared spectrum, the most prominent of which corre- 
spond to the 0(1020) and /2(1525) resonances. These and other structures in m?{K^K~) 
are now examined. 

The signal function is given by the coherent sum over resonant states that decay into 
K^K^ , plus a possible non-resonant S-wave contribution^ 

2 



S{Si2, S23, 9j/i)) = 



A=0,±1 



(4) 



where A^'^{si2, 823,6 j/^) describes the decay amplitude via an intermediate resonance 
state Ri with helicity A. Note that the J/ip has the same helicity as the intermediate 
K^K~ resonance. Each Ri has an associated amplitude strength a^' and a phase 0^' for 
each helicity state A. The amplitude for resonance i?, for each is given by 

Af{s,2,S23,eji^) = An{s23) 4^"^ T,(^^^)(^)^" f^)^" 0a(^j/^), (5) 



^The interference terms between different helicities are zero because we integrate over the angular 
variable x- 
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Figure 6: Distribution of m?{K^ K^) versus m?{J/tljK^) for candidate decays within ±20 
MeV of the mass. The horizontal bands result from the 0(1020) and /2(1525) resonances. 



where Pr is the momentum of either of the two kaons in the di-kaon rest frame, tub is 
the mass, Pb is the magnitude of the J/tp three-momentum in the rest frame, 
and F^^^^ and F^"^ are the 5° meson and Ri resonance decay form factors. The orbital 
angular momenta between the J/ip and K^K~ system is given by Lb-, and the orbital 
angular momentum in the K^K~ decay is given by Lb,'-, the latter is the same as the spin 
of the K^K~ system. Since the parent 5° has spin-0 and the J/tp is a vector, when the 
K^K~ system forms a spin-0 resonance, = 1 and Lb = 0. For K^K~ resonances 
with non-zero spin. Lb can be 0, 1 or 2 (1, 2 or 3) for Lb = 1(2) and so on. We take 
the lowest Lb the default value and consider the other possibilities in the systematic 
uncertainty. 

The Blatt-Weisskopf barrier factors F^^^^ and F^"-^ 
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are 



(6) 



ir(2) 



+ 320 + 9 
Vz^ + 3z + 9 ' 



For the B meson z = v'^Pb, where r, the hadron scale, is taken as 5.0 GeV ^; for the R 



resonance z = t'^Pbi s-iid ^ is taken as 1.5 GeV ^ 26 . In both cases Zq = t^Pq where Pq 



is the decay daughter momentum at the pole mass; for the B^ decay the J/ip momentum 
is used, while for the R resonances the kaon momentum is used. 
In the helicity formalism, the angular term, T\{9kk) is defined as 



T^{eKK) = di,{eKK), (7) 
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where d is the Wigner (i-function, J is the resonance spin, O^k is the hehcity angle of the 
in the K^K~ rest frame with respect to the K^K~ direction in the 5^ rest frame, 
and may be calculated directly from the other variables as 

[m^{J/i)K+) -m'^{J/i)K-)\m{K+K-) ,, 
cos6^i^i^ = AP P ^ • 

The J/ip hehcity dependent term 0a(6'j/^) is defined as 



Q\{0.j/^) = \Jsm^9j/^ for A = 

= yilif^ for |A| = 1. (9) 

The mass squared shape of each resonance, R is described by the function Ar{s2^)- 
In most cases this is a Breit- Wigner (BW) amplitude. When a decay channel opens close 
to the resonant mass, complications arise, since the proximity of the second threshold 
distorts the line shape of the amplitude. The /o(980) can decay to either tttt or KK. 
While the tttt channel opens at much lower masses, the K^K~ decay channel opens near 



the resonance mass. Thus, for the /o(980) we use a Flatte model 27 that takes into 
account these coupled channels. 

We describe the BW amplitude for a resonance decaying into two spin-0 particles, 
labeled as 2 and 3, as 

Ms2^) = —2 — , (10) 

where mn is the resonance mass, r(s23) is its energy-dependent width that is parametrized 
as 

rW^r.f#LV'"7^Vl. (n) 



Here Fq is the decay width when the invariant mass of the daughter combinations is equal 
to mR. 

The Flatte mass shape is parametrized as 

Ms2z) = — . r 7 -y (12) 

- S23 - imR{g^.^p^^ + gKKPKK) 

where the constants q-^t, and qkk are the /o(980) couplings to ti^ti~ and K^K~ final 
states, respectively. The p factors are given by Lorentz-invariant phase space 



P.^ = TTt/l - + ^\^ (13) 

1 / 4m?.± 1 / Amlo , 

PKK = oWi-— ^ + o\ i--:^- (14) 
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For non-resonant processes, the amphtude ^(si2, S23, is constant over the vari- 
ables S12 and S23, but has an angular dependence due to the J /ip decay. The amplitude is 
derived from Eq. ([s]), assuming that the non-resonant K^K~ contribution is an S-wave 
(i.e. Lpt = 0, Lb = I) and is uniform in phase space (i.e. Ar = 1), 



Pp. 



A{si2, S23, Oj/^) = —\ sin^ Oj/^. (15) 
4.2 Detection efficiency 

The detection efficiency is determined from a phase space simulation sample containing 
3.4 X 10^ -Bg — 7- J/iIjK^K~ events with J/ip ^ A*'*'/^"- We also use a separate sample 
of 1.3 X 10^ — )■ J/ip(f) events. The p and px distributions of the generated mesons 
are weighted to match the distributions found using J/ip(l) data. The simulation is also 
corrected by weighting for difference between the simulated kaon detection efficiencies and 
the measured ones determined by using a sample of D*~^ — )■ 7r+(i5° — > K'n'^) events. 

Next we describe the efficiency in terms of the analysis variables. Both S12 and S13 
range from 12.5 GeV^ to 24.0 GeV^, where S13 is defined below, and thus are centered at 
So = 18.25 GeV^. We model the detection efficiency using the dimensionless symmetric 
Dalitz plot observables 

X = (si2 - So) / (1 GeV^) , y = {s^, - So) / (l GeV^) , (16) 

and the angular variable 6j/^. The observables S12 and S13 are related to S23 as 

S12 + si3 + S23 = m% + mj/^ + + m|-- . (17) 

To parametrize this efficiency, we fit the cos6j/^ distributions of the J/ipK'^K~ and 
J/ip4' simulation samples in bins of m^{K^K^) with the function 

^2(S23,^J/V.) = ^^^^^^/^, (18) 

giving values of a as a function of rn?[K^ K"). The resulting distribution, shown in Fig. [7| 
is described by an exponential function 

a(s23) = exp(ai -|- a2S23), (19) 



with ai = -0.76 ± 0.18 and 02 = (-1.02 ± 0.15) GeV^^. Equation (|l8| is normalized 
with respect to cos9j/^. The efficiency in cos9j/^ depends on S23, and is observed to be 
independent of S12. Thus the detection efficiency can be expressed as 

£^(Sl2, S23, dj/i,) = y) X £:2(S23, dj/^). (20) 



After integrating over cosOj/^, Eq. (20) becomes 

"+1 



£{si2,S23,0j/^)dcosej/^ = ei{x,y), (21) 
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Figure 7: Exponential fit to the efficiency parameter a(s23)- The point near the i;^(1020) meson 
mass is determined more precisely due to the use of a large simulation sample. 



and is modeled by a symmetric fifth order polynomial function given by 

ei{x, y) = 1 + e[{x + y) + e'^ix + + e'^xy + e^{x + yf + e'^xy{x + y) 
+4(2^ + vY + (^'jxyi^x + yf + e'^x^y"^ + e'g{x + yf 
+e\^xy{x + yf + e\^xY{x + y), (22) 

where are the fit parameters. The 5° — )■ J/ijjK~^K~ phase space simulation sample 
is modeled with the polynomial function. The fitted function is shown in Fig. |8| and 
the projections of the fit are shown in Fig. |9| The efficiency is well described by the 




Figure 8: Parametrized detection efficiency as a function of m?'{K^K ) versus m?[J/ipK^). 
The z-axis scale is arbitrary. 
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Figure 9: Projections of the invariant mass squared (a) K^K^ and (b) J/ipK^ from the 
simulation used to measure the efficiency parameters. The points represent the generated event 
distributions and the curves the polynomial fit. 



parametrization. 

For the region within ±20 MeV of the 0(1020) mass, the cosOkk acceptance is used 
separately, due to the large number of signal events. Here the cos 9kk distribution shows 
a variation in efficiency, which can be parametrized using the efficiency function 



A{9_ 



KK 



1 + 

^12 COS^ GrK 
1 + e'12/3 ' 



(23) 



where the parameter is measured from a fit to the simulated J/ip(f) sample with 
ei{x,y) X A{6kk), giving = —0.099 ± 0.010, as shown in Fig. [lO 
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Figure 10: Distribution of cos 9kk for the J/')p4' simulated sample fitted with ei(x, y) x A(6kk), 
within ±20 MeV of the 0(1020) mass. 
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The mass resolution is ~ 0.7 MeV at the 0(1020) mass peak, which is added to the fit 
model by increasing the Breit-Wigner width of the 0(1020) to 4.59 MeV. 



4.3 Background composition 

The shape of the combinatorial background is modeled as 



Pb 

Cl("Sl2; "S23) \~ 

TUB 



C(Sl2,S23,6'j/^) 

where Ci(si2, S23) is parametrized as 



Co 



rrir 



823? + mlTl\ 



(l + acos'^j/^) , (24) 



Ci(si2, S23) = 1 + ci(x + y) + C2(x + yY + c^xy + cj^x + y)^ + c-^xy{x + y), (25) 



with Cj, mo, To and a as the fit parameters. The variables x and y are defined in Eq. ( |16[ ). 

shows the mass squared projections from the 5° mass sidebands with the 
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with a function of 



fit projections overlaid. The x^/iidf of the fit is 291/305. The value of a is determined 
by fitting the cos distribution of background, as shown in Fig. 
the form 1 + a cos^ Qji^-, yielding a = —0.14 ± 0.08. 
The reflection background is parametrized as 



-R(si2, S23, 9j/^) = -Ri(si2, S23) X (1 + /3 cos^ 9j/^) , (26) 

where -Ri (512,523) is modeled using the simulation; the projections of S12 and S23 are 
shown in Fig. 13 The J/ip helicity angle dependent part of the reflections is modeled 
as 1 + /3cos^ 9j/^, where the parameter (3 is obtained from a fit to the simulated cos9j/^ 
distribution, shown in Fig. 14, giving /3 = —0.19 ± 0.01. 




Figure 11: Invariant mass squared projections of (a) K^K and (b) J/ipK^ from the back- 
ground Dalitz plot of candidates in the 5^ mass sidebands. 
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Figure 12: Distribution of cos^j/^ from the background sample fit with the function 1 + 
a cos^ Q 
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Figure 13: 



Projections of the reflection background in the variables (a) m?{K^K ) and (b) 
), obtained from J/V'F*°(892) and ^° J/^K*2{U30) simulations. 



5 Final state composition 
5.1 Resonance models 

The resonances that are likely to contribute are produced from the ss system in Fig. [T| 
and thus are isoscalar (J = 0). The K^K~ system in the decay — )■ J /ipK^K~ 
can, in principle, have zero or any positive integer angular momentum. Both the P- 
parity and C-parity of K^K~ pair in a state of relative angular momentum L are given 
by (—1)^. Therefore the allowed resonances decaying to K^K~ are limited to J^^ = 
0^+, 1 , 2++, with isospin / = 0. In the kinematically accessible mass range up to 
2 GeV, resonances with J^'-' = 3 or higher are not expected and thus the subsequent 
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Figure 14: Distribution of cosOjj^ for the reflection fit with the function 1 + /3cos^ ^j/i/;- 

analysis only uses spins up to J = 2. Possible resonance candidates are listed in Table [l] 
There could also be a contribution from non-resonant events which we assume to be 
S-wave and evenly distributed over the available phase space. 



Table 1: Possible resonance candidates in the i?^ — J/ipK^K decay mode 



Spin 


Helicity 


Resonance 




Amplitude 








/o(980) 




Platte 








/o(1370), /o(1500), 


/o(irio) 


BW 


1 


0,±1 


0(1020), 0(1680) 




BW 


2 


0,±1 


/2(1270), /^(1525), 


/2(1640), 


BW 






/2(1750), /2(1950) 







To study the resonant structures of the decay 5° — )■ J/ipK^K' we use 20,425 can- 
didates with an invariant mass within ±20 MeV of the observed S° mass peak. This 
includes both signal and background, with 94% signal purity. We begin our analysis con- 
sidering only the resonance components 0(1020), /2(1525) and a non-resonant component, 
established in our earlier measurement [lOj , and add resonances until no others are found 
with more than two standard deviation statistical significance (2cr). The significance is 
estimated from the fit fraction divided by its statistical uncertainty. Our best fit model 
includes a non-resonant component and 8 resonance states: 0(1020), /o(980), /o(1370), 
/^(1525), /2(1640) (|A| = 1), 0(1680) (|A| = 1)[] /2(1750), and /2(1950). Most of the res- 
onances considered here are well established except for the modes /2(1640), /2(1750), and 
/2(1950). Although the existence of /2(1640) is not confirmed yet [9], the right shoulder 



"The /2(1640) (A = 0) and (/)(1680) (A — 0) components have less than two standard deviation 
significance when added separately to the fit, and therefore are not included in the best fit model. 
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of /2(1525) fits better when we add this state. The presence of multiple broad overlapping 
resonances in this region may indicate a failure of the isobar model used in this analysis, 
but with the present data sample alternative descriptions are not feasible. Indeed, the 
situation is not clear for the resonance states in the vicinity of 1750 MeV. The PDG 
lists a spin-0 resonance, /o(1710), around 1.72 GeV of K^K" invariant mass |9j. The 
Belle collaboration observed a resonance in the vicinity of 1.75 GeV with J^^ = (even)^"'" 



in their study of 77 — )• K^K~ 28 , but could not establish its spin. A state of mass 
1767±14 MeV was seen by the L3 collaboration decaying into K^K^ with J = 2 [29]. 
We find that our data are better fit including the /2(1750) mode. If we substitute either 
the /o(1710) or /o(1750) resonance the fit is worsened, as the — ln£ increase by 59 and 7 
units, respectively. 

In the same analysis of 77 — )■ K^K~, Belle also observed the /2(1950) [28] resonance. 
We include this state in our best fit model. Furthermore, we do not expect significant con- 
tributions from the /2(1270) and /o(1500) resonances, since the PDG branching fractions 
are much larger in the vr+vr" final state than in K^K~ j9] and we did not see significant 
contributions from these two resonances in the — >■ J/ip7r~^n~ final state |8]. Therefore, 
these two resonances are not considered in the best fit model. However, we add these 
states, in turn, to the best fit model in order to test for their possible presence. 

The masses and widths of the BW resonances are listed in Table [H When used in the 
fit they are fixed to the central values, except for the /2(1525), whose mass and width are 
allowed to vary. 



Table 2: Breit-Wigner resonance parameters. 
Resonance Mass (MeV) Width (MeV) Source 



0(1020) 
/2(1270) 
/o(1370) 
/o(1500) 
/^(1525) 
/2(1640) 
0(1680) 

/o(irio) 

/2(1750) 
/2(1950) 



1019.46 ±0.02 
1275±1 
1475±6 
1505±6 
1525±5 
1639±6 
1680±20 
1720±6 
1737±9 
1980±14 



4.26±0.04 
185±3 
113±11 
109±7 
73±6 
99±60 
150±50 
135±8 
151±33 
297±13 



PDG 
PDG 
LHCbT8 
PDG 
PDG 
PDG 
PDG 
PDG 
Belle 
Belle 
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The /o(980) is described by a Flatte resonance shape, see Eq. (12). The parameters 
describing the function are the mass, and the couplings Qt^t^ and Qkk, which are fixed in the 
fit from the previous analysis of — )■ J/ipTc^iT^ |8|. The parameters are rriQ = 939.9±6.3 
MeV, g^TT = 199 ± 30 MeV and gKK/g-n-K = 3.0 ± 0.3. All background and efficiency 
parameters are fixed in the fit. 

To determine the complex amplitudes in a specific model, the data are fitted maxi- 
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mizing the unbinned likelihood given as 



N 



125 -523' ^J/ip)y 



(27) 



where is the total number of candidates, and F is the total PDF defined in Eq. ([T|. 
The PDF normalization is accomplished by first normalizing the J/ ip helicity dependent 
part by analytical integration, and then for the mass dependent part using numerical 
integration over 400x800 bins. 



The fit determines the relative values of the amplitude strengths. 



defined in Eq. (4). We choose to fix a, 



0(1020) 



-.Ri 



and phases. 



1. As only relative phases are 



physically meaningful, one phase in each helicity grouping must be fixed. In addition, 
because J/ipK^K' is a self-charge-conjugate mode and does not determine the initial 
B fiavor, the signal function is an average of and 5°. If we consider no K^K~ 
partial-waves of a higher order than D-wave, then we can express the differential decay 
rate (dV /drriKK dcosOKK dcosOj/^) derived from Eq. ^ in terms of S-, P-, and D-waves 
including helicity and ±1 components. The differential decay rates for B^ and B^, 
respectively are 

df 



drriKK d cos 9kk d cos 9j/^ 



A'so^"'^ + ^Poe^'"° cos 9kk + A^/"'^^^ ( ^ cos^ 9kk 



sm^9 



J/V 



+ 



and 



^p±ie""±^ \l I sin 9kk + ^Li^e^'^^ii J ^ sin 9kk cos 9kk 



dr 



1 + cos^ 9 



(28) 



drriKK d cos 9kk d cos 9 j/^ 



A's/^'o - ^Po^^""" cos^KK + ^Loe ( ^ cos^ 9kk 



sin^ 9j/^ 



+ 



A'p^^e^p±^ J I sin 9kk - A^^/'''''^^ J ^ sin 9kk cos 9kk 



1 + cos^ 9 



(29) 



where Al^ and (pl^ are the sum of amplitudes and reference phases, for the spin-fc reso- 
nance group, respectively. The decay rate for B^ is similar to that of 5°, except 9k+k- 
and 9ji^ are now changed to tt — 9k+k- and vr — 9ji^ respectively, as a result of using K~ 
and to define the helicity angles and hence the signs change in front of the A^p^ and 
A]^^^ terms. 



Summing Eqs. (28) and (29) results in cancellation of the interference involving A = 



terms for spin-1, and the A = ±1 terms for spin-2, as they appear with opposite signs for 
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and decays. Therefore we have to fix one phase in the spin-1 (A = 0) group 



and one in the spin-2 (A = ±1) group 



The other phases in each corresponding 



group are determined relative to that of the fixed resonance. 



5.2 Fit results 

The goodness of fit is calculated from 3D partitions of Si2, S23 and cos9j/^. We use the 



Poisson hkelihood x 30 defined as 



N 



X 



i=l 



Hi + njln 



(30) 



where is the number of candidates in the three dimensional bin i and Xj is the expected 
number of candidates in that bin according to the fitted likelihood function. An adaptive 
binning algorithm is used, requiring a minimum of 25 entries in each bin. The associated 
number of degrees of freedom (ndf ) is N — k — 1, where k is the number of free parameters 
in the likelihood function. The x^/ndf and the negative of the logarithm of the likelihood, 
— Ini2, of the fits are given in Table [3j Starting values of parameters are varied in order 
to ensure that global likelihood minimums are found rather than local minimums. 



Table 3: /jidi and — ln£ of different resonance models. 



Resonance model 


-ln£ 


X^/ndf 


Best fit 
Best fit + /2(1270) 
Best fit + /o(1500) 


29,275 
29,273 
29,274 


649/545 = 1.1908 
644/541 = 1.1911 
647/543 = 1.1915 



Attempts to add one more resonance such as /2(1270) and /o(1500) improve the — ln£ 
marginally, but the x^/ndi are worse than the best fit model. We retain only those 
resonances that are more than 2a significant, except for the /2(1750) where we allow the 
|A| = 1 component, since the A = component is significant. For models with one more 



resonance, the additional components never have more than 2a significance. Figure 15 
shows the projection of im?{K^K~) for the best fit model, the Tm?{J /ipK^) and cosOj/^ 



projections are displayed in Fig. 16 The projection of the K'^K' invariant mass spectrum 
is shown in Fig. [TT} 

While a complete description of the 5° — )■ J/ipK'^K~ decay is given in terms of the 
fitted amplitudes and phases, knowledge of the contribution of each component can be 
summarized by defining a fit fraction, J^^. To determine J^^ we integrate the squared 
amplitude of R over the Dalitz plot. The yield is then normalized by integrating the entire 
signal function over the same area. Specifically, 



-cR 
■J X 



dsi2 ds23 dcosOj/^ 



J S{si2,S23,0j/^) dsi2 ds23 dcosO 



J/4, 



(31) 
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Figure 15: Dalitz plot fit projection of m^{K^K'') using a logaritlimic scale. The points with 
error bars are data, the (black) dotted curve shows the combinatorial background, the (red) 
dashed curve indicates the reflection from the misidentified — )• J/7pK~TT~^ decays, and the 
(blue) solid line represents the total. 




Figure 16: Dalitz plot fit projections of (a) m?{J/ipK~^) and (b) cosOj/^. The points with error 
bars are data, the (black) dotted curve shows the combinatorial background, the (red) dashed 
curve indicates the reflection from the misidentified — ?• J/tpK~iT~^ decays, and the (blue) 
solid line represents the total fit results. 
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Figure 17: Dalitz fit projection of m{K^K~). The points represent the data, the dotted 
(blaclc) curve shows the combinatorial background, and the dashed (red) curve indicates the 
reflection from misidentified — )■ J/ijjK~Tr~^ decays. The largest three resonances 4>{W20), 
/2(1525) and /o(980) are shown by magenta, brown and green long-dashed curves, respectively; 
all other resonances are shown by thin black curves. The dashed (cyan) curve is the non-resonant 
contribution. The dot-dashed (black) curve is the contribution from the interferences, and the 
solid (blue) curve represents the total fit result. 



Note that the sum of the fit fractions is not necessarily unity due to the potential presence 
of interference between two resonances. Interference term fractions are given by 



-RR' ^ I j a'^ e^'^'^ ^A^{su,S23,9j/^)A^'*{si2,S23,0j/^)dsi2 ds23 dcosOj/^ 



J-f « = 7^e 



/ S{si2,S23,0j/^) dsi2 ds23 dcos9j/^ 



(32) 



and 



E E-^a^+E-^a""' =1. (33) 

A \ R RR' ) 

If the Dalitz plot has more destructive interference than constructive interference, the 
sum of the fit fractions will be greater than unity. Conversely, the sum will be less than 
one if the Dalitz plot exhibits constructive interference. Note that interference between 
different spin-J states vanishes because the d^^^^ angular functions in A^ are orthogonal. 

The determination of the statistical uncertainties of the fit fractions is difficult because 
they depend on the statistical uncertainty of every fitted magnitude and phase. Therefore 
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we determine the uncertainties from simulated experiments. We perform 500 experiments: 
eacli sample is generated according to tlie model PDF, input parameters are taken from 
the fit to the data. The correlations of fitted parameters are also taken into account. 
For each experiment the fit fractions are calculated. The distributions of the obtained 
fit fractions are described by Gaussian functions. The r.m.s. widths of the Gaussian 
functions are taken as the statistical uncertainties on the corresponding parameters. The 
fit fractions and phases of the contributing components are given in Table [4], while the fit 
fractions of the interference terms are quoted in Table [5] 



Table 4: Fit fractions (%) and phases of contributing components. For P- and D-waves 
A represents the helicity. 



Component Fit fraction (%) Phase (degree) 



0(1020), 


A = 





32.1 ±0.5 ±0.8 


(fixed) 


0(1020), 


|A| = 


= 1 


34.6 ±0.5 ±1.3 


(fixed) 


/o(980) 






12.0 ± i.stli 


-294 ± 8 ± 25 


/o(1370) 






1.2 ±0.3^?:^ 


-81 ±8±8 


/^(1525) 


A = 





9.9 ±0.71^-4 


(fixed) 


/^(1525) 


|A| 


= 1 


5.1 ±0.9l}i 


(fixed) 


/2(1640) 


|A| 


= 1 


1.5±0.7^°;9^ 


165 ± 27tll 


0(1680), 


|A| = 


= 1 


3.4±0.3l^:^ 


106 ± 141^60 


/2(1750) 


A = 





2.6±0.5lJ:^ 


238 ± 8 ± 9 


/2(1750) 


|A| 


= 1 


1.8±1.0l^j 


45 ± 30t]l 


/2(1950) 


A = 





0.4±0.2+°i 


46 ± 17^^^° 


/2(1950) 


|A| 


= 1 


1.7 ±0.5^^:^ 


-53 ± 26^1^° 


Non-resonant 




6.0 ±1.6^1^ 


-39 ± 6 ± 23 


Total 






112.3 





Table [6] shows a comparison of the fit fractions when a different parametrization is used 
for the /o(980) resonance shape. The BES /o(980) functional form is the same as ours with 
the parameters tuq = 965 ±10 MeV, g^^ = 165 ±18 MeV and gKK/gnn = 4.21 ±0.33 |31]. 
The BaBar collaboration assumes that the non-resonant S-wave is small and consistent 
with zero |32j. The BaBar functional form is different and parametrized as 



1 



Ms23) = -2 : — F^— > (34) 



with pkk = 2Pr/^/s^. The parameters are rriR = 922 ± 3 MeV and = 240 ± 80 MeV, 



taken from BaBar's Dalitz plot analysis of Df — t- K^K~Tr~^ 33 . The /o(980) fraction is 
smaller in the BaBar parametrization, while the total S-wave fraction is consistent in the 
three different parameterizations. In all cases, the dominant component is the 0(1020), 
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Table 6: Comparison of the fit fractions (%) with the LHCb, BES and BaBar /o(980) 
parameterizations described in the text. For P- and D- waves, A represents the hehcity. 



Component 


LHCb 


BES 


BaBar 


0(1020), A = 


32.1 ± 0.5 


32.1±0.5 


32.0 ± 0.5 


0(1020), |A = 1 


34.6 ± 0.5 


34.6 ± 0.5 


34.5 ± 0.5 


/o(980) 


12.0 ± 1.8 


9.2 ± 1.4 


4.8 ± 1.0 


/o(1370) 


1.2 ± 0.3 


1.2 ± 0.3 


1.3 ± 0.3 


/^(1525), A = 


9.9 ± 0.7 


9.8 ± 0.7 


9.5 ± 0.7 


/2(1525), |A| = 1 


5.1 ± 0.9 


5.1 ± 0.9 


4.9 ± 0.9 


/2(1640), |A| = 1 


1.5 ± 0.7 


1.5 ± 0.7 


1.5 ± 0.7 


J, / -1 (2QC\\ \\\ 1 

0(^iooUj, A = 1 


3.4 ± U.3 


3.4 ± U.3 


O A _i_ n Q 
3.4 ± U.3 


/2(1750), A = 


2.6 ±0.5 


2.5 ±0.5 


2.2 ±0.5 


/2(1750), |A| = 1 


1.8 ± 1.0 


1.8 ± 1.0 


1.9 ± 1.0 


/2(1950), A = 


0.4 ±0.2 


0.4 ±0.2 


0.4 ±0.2 


/2(1950), |A| = 1 


1.7±0.5 


1.8 ±0.5 


1.8 ±0.5 


Non-resonant S-wave 


6.0 ± 1.4 


4.7 ± 1.2 


8.6 ± 1.7 


Interference between S-waves 


-5.5 


-1.7 


-1.1 


Total S-wave 


13.7 


13.4 


13.6 


-ln£ 


29,275 


29,275 


29,281 


xVndf 


649/545 


653/545 


646/545 



the second largest contribution is the /2(1525), and the third the /o(980) resonance. There 
are also significant contributions from the /o(1370), /2(1640), 0(1680), /2(1750), /2(1950) 
resonances, and non-resonant final states. The amount of /o(980) is strongly parametriza- 
tion dependent, so we treat these three models separately and do not assign any systematic 
uncertainty based on the use of these different /o(980) shapes. Therefore we refrain from 
quoting a branching fraction measurement for the decay 5° — )■ J/ipfoldSO). 

The determination of the parameters of the /2(1525) resonance are not dependent on 
the /o(980) parametrization. The parameters of the /2(1525) are determined to be: 

m/^(i525) = 1522.2 ± 2.8^11^ MeV, 
r/^(i525) = 84 ±6t^°MeV. 

Whenever two or more uncertainties are quoted, the first is the statistical and the second 
systematic. The latter will be discussed in Section [5l6| These values are the most accurate 
determinations of the /2(1525) resonant parameters [9]. Note that our determination of 
the mass has the same uncertainty as the current PDG average. 
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5.3 K^K S-wave in the 0(1020) mass region 



It was claimed by Stone and Zhang [12] that in the decay of 5° — )■ J/ip4>-, the K^K~ 
system can have S-wave contributions under the 0(1020) peak of order 7% of the total 
yield. In order to investigate this possibility we calculate the S-wave fractions as given 
by the fit in 4 MeV mass intervals between 990 < m{K^K~) < 1050 MeV. The resulting 
behavior is shown in Fig. 18 Here we show the result from our preferred model and 
also from the alternative /o(980) parameterizations discussed above. The observation of 
significant S-wave fractions in this region means that this contribution must be taken into 
account when measuring CP violation in the (f) mass region. The total S-wave fraction 
as a function of the mass interval around the (j) mass is also shown in Fig. 19 Using 
a time dependent analysis of 5° J/^/'0(lO2O), LHCb reported (2.2 ± 1.2 ± 0.07)% § 
of S-wave within ±12 MeV of the 0(1020) mass peak. We measure the S-wave fraction 
within the same mass window as a consistent, and more precise (1.1 ± 0.11q^)%. CDF 
measured the S-wave fraction as (0.8 ± 0.2)% for m{K~^K~) within about ±9.5 MeV of 
the mass [6], while ATLAS quotes (2±2)% for an 11 MeV interval These results 
are consistent with ours. The DO collaboration, however, claimed a (14.7±3.5)% S-wave 
fraction within approximately ±10 MeV of the meson mass js], in disagreement with 
all of the other results. 
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Figure 18: S-wave fraction as a function of m{K^K^) starting from 990 MeV up to 1050 MeV 
in 4 MeV mass intervals. The squares (blue), triangles (red), and circles (green) represent the 
LHCb, BES and BaBar parameterizations of /o(980), respectively. The experimental statistical 
uncertainties are only shown for the LHCb model; they are almost identical for the other cases. 
The experimental mass resolution is not unfolded. 
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Figure 19: S-wave fractions in different m{K~^ K~) intervals centered on the (p meson mass. 
The squares (blue), triangles (red), and circles (green) represent the LHCb, BES and BaBar 
parameterizations of /o(980), respectively. The experimental statistical uncertainties are only 
shown for the LHCb model; they are almost identical for the other cases. The experimental 
mass resolution is not unfolded. 



5.4 Helicity angle distributions 

The decay angular distributions or the helicity angle distributions are already included in 
the signal model via Eqs. (jsj) and ([9|. In order to test the fit model we examine the cos Oj/^ 
and cosOkk distributions in two different K^K~ mass regions: one is the 0(1020) region 
defined within ±12 MeV of the 0(1020) mass peak and the other is defined within one full 
width of the /2(1525) mass. The background-subtracted efficiency-corrected distributions 
are shown in Figs. |20] and 21 The distributions are in good agreement with the fit model. 



5.5 Angular moments 

The angular moment distributions provide an additional way of visualizing the effects 
of different resonances and their interferences, similar to a partial wave analysis. This 
technique has been used in previous studies [8,34 . 

We define the angular moments (Yl^) as the efficiency-corrected and background- 
subtracted K^K~ invariant mass distributions, weighted by orthogonal and normalized 
spherical harmonic functions YJ°(cos ^^x), 

(Yi^) = J dT{mKK, cos 9kk)Yi° {cos 9KK)dcos9KK- (35) 

If we assume that no K'^K' partial- waves of a higher order than D-wave contribute, 
then we can express the differential decay rate, derived from Eq. (|4]) in terms of S-, P-, 
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Figure 20: Background-subtracted efficiency-corrected cos^j/^ helicity distributions: (a) in 
(/>(1020) mass region (x^/ndf = 54.4/40), (b) in /2(1525) mass region (x^/ndf = 34.4/40). 




Figure 21: Background-subtracted efficiency-corrected cosOkk helicity distributions: (a) in 
(/)(1020) mass region (xV^df = 57.4/40), (b) in /2(1525) mass region {x^/ndf= 43.4/40). The 
distributions are compatible with expectations for spin-1 and spin-2, respectively. 



and D-waves including helicity and ±1 components as 
dT 



drriKK d cos 9kk 



271 \SoY^\cos9kk) + Poe**^«Fi°(cos^i^K) + Voe"^°'^Y^icos9KK 



+ 27T 




3 /IS 

sin Okk + V±ie'^''±^ \ — sin Ork cos 6kk 
Sir V oTT 



(36) 



where S\, Vx, V\ and are real- valued functions of rriKK, and we have factored out 
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the S-wave phase. 



We can then calculate the angular moments 



25oPoCOS<l>p„ + 4=n^^0COs(<l>p„ - <^Do) 

V5 



4vr(F3°) 



2^5 



(37) 



Vl + 2S^V^ cos$B„ + 



6\/-Po^^oCos($p„-$z)o, 
6^2 4 . 



1 




7 

6 



7 



Piir'il COs($p±, -'^'Dii; 



The angular moments for / > 4 vanish. Figures 22 and 23 show the distributions of the 
angular moments for the fit model around ±30 MeV of the 0(1020) mass peak and above 
the 0(1020), respectively. In general the interpretation of these moments is that (Xq) is 
the efficiency-corrected and background-subtracted event distribution, (Y^) the sum of the 
interference between S-wave and P-wave, and P-wave and D-wave amplitudes, (1^2*^) the 
sum of P-wave, D-wave and the interference of S-wave and D-wave amplitudes, (1^3°) the 
interferenc e be tween P-wave and D-wave amplitudes, and (1^4°) the D-wave. As discussed 
the average of and 5° cancels the interference terms that involve Vq 



5.1 



in Section 

and This causes the angular moments (Yi) and i^^) to be zero when averaging over 
-B^ and -B^ decays. We observe that the fit results well describe the moment distributions, 
except for the (Yf) and (y^) values below 1.2 GeV. This may be the result of statistical 
fiuctuations or imperfect modeling. 



5.6 Systematic uncertainties 

The sources of the systematic uncertainties on the results of the Dalitz plot analysis are 
summarized in Table [7} The uncertainties due to the background parametrization are es- 
timated by comparing the results from the best fit model with those when the background 
shape parameters are obtained from a fit to the lower sideband region only. The uncer- 
tainties in the efficiency are estimated by comparing the fit results when the efficiency 
parameters are changed by their statistical uncertainties and are added in quadrature. 
The effect on the fit fractions of changing the efficiency function is evaluated using a 
similar method to that used previously Briefiy, we change the efficiency model by 
increasing the minimum IP requirement from 9 to 12.5 on both of the kaon candidates. 
This has the effect of increasing the of the fit to the angular distributions of i?^ — t- J/il) cf) 
data by 1 unit. The new efficiency function is then applied to the data with the original 
minimum IP selection of 9, the likelihood is re-evaluated and the uncertainties are 
estimated by comparing the results with the best fit model. The largest variations among 
these two efficiency categories are included in the uncertainty. 
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Figure 22: Dependence of the spherical harmonic moments of cosOkk as a function of the 
K^K~ mass around the (/)(1020) mass peak after efficiency corrections and background sub- 
traction. The points with error bars are the data and the solid curves are derived from the fit 
model. 



We estimate additional uncertainties by comparing the results when one more res- 
onance is added to the best fit model. The uncertainties due to the line shape of the 
contributing resonances with fixed mass and width parameters are estimated by varying 
them individually in the fit according to their combined statistical and systematic uncer- 
tainties added in quadrature. We compare the results with the best fit and add them in 
quadrature to estimate the uncertainties due to the hne shape. 

Another source of systematic uncertainty is the value we choose for Lb-, the orbital 
angular momentum in the decay. If Lr equals zero then Lb equals zero. If, however, 
Lb, is 1 then Lb can either be or 1, and if Lb is 2, Lb can be 1, 2 or 3. For our 
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Figure 23: Dependence of the spherical harmonic moments of cos 9kk as a function of the 
K^K~ mass above 1050 MeV, after efficiency corrections and background subtraction. The 
points with error bars are the data and the sohd curves are derived from the fit model. 



best fit we don't allow multiple values for Lb, but choose the lowest allowed value. To 
estimate the systematic uncertainties due to the choice of Lb, we repeat the fit changing 
the default value of Lb, in turn, to each higher allowed value and compare the fit results 
with the best fit. The differences are grouped into the fit model category, and we assign 
the largest variations as the systematic uncertainties. These later two categories often 
give in asymmetric uncertainties. 
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Table 7: Absolute systematic uncertainties on the fit results. 



Item 


Efficiency 


Background 


Fit model 


Total 


^/^(i525) (MeV) 


1.2 


0.4 


+5.2 
-1.5 


+5.3 
-2.0 


r/^(i525) (MeV) 


4.7 


0.5 


+8.6 
-1.8 


+9.8 
-5.0 


Fit fractions (%) 


0(1020) A = 


0.8 





+0.06 

— U.U4 


±0.8 


0(1020) lAl = 1 


1.3 





+0.4 
— U.l 


±1.3 


fo(980) 


1.7 


0.4 


+2.1 
—1.7 


+2.8 
—2.5 


/o(1370) 


0.3 


0.02 


+0.2 

1 9 


+0.3 

— 1 .0 


/^(1525) A = 


1.5 


0.2 


+ 1.8 

—0 5 


+2.4 
— 1.6 


/^(1525) |A| = 1 


1.1 


0.4 


+ 1.3 
—0.8 


+1.8 
— 1.4 


/2(1640) |A| = 1 


0.1 


0.1 


+0.7 
—0.9 


+0.7 
—0.9 


0(1680) |A| = 1 


0.3 


0.1 


+4.4 
—0.1 


+4.4 
—0.3 


/2(1750) A = 


0.3 


0.1 


+1.0 
—0.5 


+1.0 
—0.6 


/2(1750) |A| = 1 


1.6 


0.3 


+1.5 
—2.9 


+2.2 
—3.3 


/2(1950) A = 


0.1 


0.04 


+0.2 


+0.2 


/2(1950) |A| = 1 


2.1 


0.6 


+1.1 

—3.1 


+2.5 
—3.8 


Non-resonant 


1.7 


0.4 


+0.9 
—1.5 


+2.0 
—2.2 


S-wave within ±12 McV 


0.02 


0.02 


+0.2 
—0.1 


+0.2 
—0.1 


of 0(1020) peak (%) 












Phases (dej 


^rees) 






/o(980) 


19 


8 


+14 
- 8 


+25 
-22 


/o(1370) 


6 


1 


+6 

-5 


±8 


/2(1640) |A| = 1 


7 


11 


+ 
-42 


+13 
-44 


0(1680) |A| = 1 


1 


1 


+260 
-210 


+260 
-210 


/2(1750) A = 


2 


2 


+9 

-9 


±9 


/2(1750) |A| = 1 


12 


10 


+ 2 
-70 


+16 
-70 


/2(1950) A = 


13 


1 


+108 
- 14 


+110 
- 20 


/2(1950) |A| = 1 


77 


26 


+130 
- 1 


+150 
- 80 


Non-resonant 


18 


7 


+12 
- 5 


+23 
-20 



29 



6 Absolute branching fractions 



Branching fractions are measured from ratios of the decay rates of interest normahzed to 
the well established decay mode B~ J/ipK'. This decay mode, in addition to having 
a well measured branching fraction, has the advantage of having two muons in the final 
state and hence the same triggers as the decay. However, we require knowledge of 
the B^/B~ production ratio. For this we assume isospin invariance and use the B^/B^ 
production ratio fs/fd = 0.256 ± 0.020, given in Ref. [35]. The branching fractions are 
calculated using 

^ J/^X) = X BiB- ^ J/i^K-) X -i-, (38) 

Nb-I^b- J slid 

where X indicates a specific K^K~ state, represents the yield of the decay of interest, 
and e corresponds to the overall efficiency. We form an average of B{B~ — )■ J/ipK') = 



(10. 18 ±0.42) X 10 ^ using the recent Belle |36] and BaBar 37 measurements, corrected 
to take into account different rates of B^B~ and B^B^ pair production from T{4S) using 
^g^ = 1.055 ±0.025 jgj. 

The detection efficiency is obtained from simulation and is a product of the geometrical 
acceptance of the detector, the combined reconstruction and selection efficiency and the 
trigger efficiency. The efficiency also includes the efficiency of the Dalitz plot model for 
the case of B^ — )■ J/tpK^K" , where the best fit model is used. The detection efficiencies 
and their various correction factors are given in Table |8} To ensure that the p and px 
distributions of the generated B meson are correct we weight the B^ simulations using 

— >■ J/'?/'0(lO2O) data and the B~ simulations using B~ — )• J/ipK~ data. Since the 
control channel has a different number of charged tracks than the decay channel, we weight 
the simulations with the tracking efficiency ratio by comparing the data and simulations 
in bins of the track's p and px- we further weight the B^ — )■ J/ipK~^K~ simulation, using 
the PDG value of 5° lifetime, (1.497 ± 0.015) x lO'^^ g ^^^^^ 



Table 8: Detector efficiencies determined from simulation and the correction factors. 



Item 


J/ijK+K- 


Jli)K- 


Detection efficiency (%) 


1.061 ±0.004 


2.978 ±0.011 


Correction factors 


Tracking efficiency 


0.999 ±0.010 


1.003 ±0.010 


PID 


0.819 ±0.008 


0.974 ±0.005 


p and pt 


1.077 ±0.005 


1.053 ±0.005 


B^ lifetime 


0.993 ±0.015 




Total efficiency (%) 


0.887 ±0.004 ±0.018 


3.065 ±0.012 ±0.038 
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The resulting branching fractions are 



B{Bl^ J/ijK+K~) = (7.70 ± 0.08 ± 0.39 ± 0.60) X 10-^ 
J/^0(1O2O)) = (10.50 ± 0.13 ± 0.64 ± 0.82) X 10-^ 
i3(5° ^ 7/^^/2(1525)) = (2.61 ± 0.20l°;^2 ± o.20) x 10-^ 

where the branching fractions i3(0(lO2O) ^ K+R-) = (48.9 ± 0.5)% and i3(/^(1525) ^ 
K^K^) = (44.4 ± 1.1)% are used [9]. Here the first uncertainty in each case is statistical, 
the second is systematic and the third reflects the uncertainty due to fs/fd- Note that 
these are the time-integrated branching fractions. Results on the polarization fractions 
of — )■ J/ip (f){1020) from a time-dependent analysis will be forthcoming in a separate 
pubhcation f38|. The ratio of -8(5° ^ J /iljf'^{lh2h)) /B{Bl J/^/'0(lO2O)) is consistent 
with our previous result |10|, DO [Tl|, and the Belle result [39]. The current PDG value of 
i3(5° ^ J/V^0(1O2O)) = (L4 ± 0.5) X 10"^ is dominated by the CDF measurement |40]. 
Our measured value is in good agreement with this measurement and also the most recent 
yet unpublished values measured by CDF |41| and Belle 139]. The Belle collaboration 
has also recently reported the branching fraction of i3(-B° — )■ J/ipK^K') |39j, where 
5° ^ J/tlj(j){1020) is excluded. 

The systematic uncertainty on the branching fraction has several contributions listed 
in Table [9} Since the branching fractions are measured with respect to B^ — ?■ J/ijjK~ 
which has a different number of charged tracks than the decays of interest, a 1% sys- 
tematic uncertainty is assigned due to differences in the tracking performance between 



Table 9: Relative systematic uncertainties on branching fractions (%). 



Item 


J/tpK+K- 


J/?/'</)(1020) 


J/V'/^(1525) 


Tracking efficiency 


1.0 


1.0 


1.0 


Material and physical effects 


2.0 


2.0 


2.0 


PID efficiency 


1.0 


1.0 


1.0 


B^ p and distributions 


0.5 


0.5 


0.5 


B~ p and px distributions 


0.5 


0.5 


0.5 


lifetime 


1.5 


1.5 


1.5 


Efficiency function 


0.02 


0.02 


0.02 


B(0(1O2O) ^ K+K-) 




1.0 




fi(/^(1525) -> 






2.5 


B{B- J/i)K-) 


4.1 


4.1 


4.1 


Contributions 


from Dalitz analysis 




Efficiency 




3.1 


12.3 


Background 







1.3 


Fit model 




+0.7 
-0.2 


+14.6 
-11.2 


Sum in quadrature of items above 


5.0 


+6.1 
-6.0 


+20.0 
-17.7 


fs/ fd 


7.8 


7.8 


7.8 
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data and simulation. Another 2% uncertainty is assigned for the additional kaon which 
is due to decay in flight, large multiple scatterings and hadronic interactions along the 
track. Using the PDG value for the lifetime [o] as input gives rise to an additional 
1.5% systematic uncertainty. Small uncertainties are introduced if the simulation does 
not have the correct B meson kinematic distributions. We are relatively insensitive to 
any of these differences in the B meson p and distributions since we are measuring 
the relative rates. By varying the p and p^ distributions we see at most a change of 
0.5%. There is a 1% systematic uncertainty assigned for the relative particle identiflca- 
tion efficiencies. An uncertainty of 0.02% is included due to the change of the efficiency 
function Eq. ( [20) ). Three additional uncertainties are considered in the branching frac- 
tions of J/ip(f){1020)) and -)■ J/V'/2(1525)) as these are measured from the 
fit fractions of the Dalitz plot analysis. The total systematic uncertainty is obtained by 
adding each source of systematic uncertainty in quadrature as they are uncorrelated. 



7 Conclusions 

We have determined the flnal state composition of the B^ — ?■ J/ipK^K' decay channel 
using a modifled Dalitz plot analysis where we include the decay angle of the J/ip. The 
largest contribution is the 0(1020) resonance, along with other S-, P- and D-wave K^K~ 
states, and a non-resonant K^K^ contribution. All of the components are listed in 
Table |4j The mass and width of the /2(1525) resonance are measured as 

m/^(i525) = 1522.2 ± 2.8111; MeV, 
r/^(i525) = 84 ±6li°MeV. 

We also observe a signiflcant S-wave component that is present over the entire K^K~ 
mass region. Within ±12 MeV of the 0(1020) mass it is (1.1 ± Qlt^DVo of the yield, and 



can affect precision CP violation measurements 12 . Finally we determine the absolute 
branching fractions 

B{Bl-^ J/i)K+K-) = (7.70 ± 0.08 ± 0.39 ± 0.60) X 10-^ 
i3(5° ^ J/^0(1O2O)) = (10.50 ± 0.13 ± 0.64 ± 0.82) X 10-^ 
i3(5° ^ J/^/^(1525)) = (2.61 ± 0.20l[!i2 ± 0.20) x 10-^ 

where the flrst uncertainty in each case is statistical, the second is systematic and the third 
due to fs/fd- These results provide a good understanding of the J/il)K'^K~ flnal state in 
-B° decays over the entire kinematically allowed region. The J/-0/2(1525) results supersede 
those of Ref. 10 . This decay mode offers the opportunity for additional measurements 
of CP violation |42 . 
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